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Applications of perturbative QCD to deeply virtual Compton scattering and hard exclusive me- 
son electroproduction processes require a generalization of usual parton distributions for the case 
when long-distance information is accumulated in nonforward matrix elements {p' \ 0{Q,z) \ p) of 
quark and gluon light-cone operators. In our previous papers we used two types of nonperturba- 
tive functions parametrizing such matrix elements: double distributions F[x, y; t) and nonforward 
distribution functions T(^{X; t). Here we discuss in more detail the double distributions (DDs) and 
evolution equations which they satisfy. We propose simple models for F{x, y;t = 0) DDs with 
correct spectral and symmetry properties which also satisfy the reduction relations connecting 
them to the usual parton densities f{x). In this way, we obtain self-consistent models for the 
^-dependence of nonforward distributions. We show that, for small one can easily obtain non- 
forward distributions (in the X > region) from the parton densities: J-q{X; t = 0) « f{X — (/2). 
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I. INTRODUCTION 



Applications of perturbative QCD to deeply virtual Compton scattering and hard exclusive electroproduction 
processes require a generalization of usual parton distributions for the case when long-distance information is 
accumulated in nonforward matrix elements {p — r\ 0(0, z) | p) | ^2^q of quark and gluon light-cone operators. As 
argued in Refs. ||^,^,^, such matrix elements can be parametrized by two basic types of nonperturbative functions. 
With z taken in the lightcone "minus" direction, the double distributions (DDs) F{x,y;t) specify the light-cone 
"plus" fractions xp~^ and yr~^ of the initial hadron momentum p and the momentum transfer r carried by the initial 
parton. Though z is an integration variable, only one direction on the lightcone (specified by external momenta) 
is important for the lightcone-dominated processes. In other words, only the lightcone plus direction of the hadron 
p and r momenta are essential for such processes. By definition, the DDs F{x,y;t) do not depend on the r^/p'^ 
ratio. On the other hand, treating the proportionality coefficient as an independent parameter: r+ = Cp^i one can 
introduce an alternative description in terms of the nonforward parton distributions J-(^(X; t) with X = x + yQ being 
the total fraction of the initial hadron momentum taken by the initial parton. The shape of the functions J^i^{X\ t) 
explicitly depends on the parameter Q characterizing the skewedness of the relevant nonforward matrix element. 
This parametrization of nonforward matrix elements by Tq{X] t) is similar to that proposed originally by X. Ji 
who introduced off-forward parton distributions (OFPDs) H{x,£,;t). The latter are close to functions considered 
earlier in Ref. |8j. The functions H{x,^;t) have a simple relation to nonforward distributions (NFPDs) J^(^{X;t), 
while the non-diagonal distributions F{xi, X2) discussed by Collins, Frankfurt and Strikman 0| essentially coincide 
with TQ=xi-x2ixilt = 0) (see Ref. for details). The basic distinction between our approach and those of Refs. 

is that we treat the double distributions F{x,y;t) as the primary objects of the QCD analysis producing 
the nonforward distributions TQ{X;t) (and other types of distributions) after an appropriate integration. 

The formalism of double distributions provides a rather effective tool for studyiiig some general (e.g., spectral) 
properties of NFPDs and it allows to find analytic solutions of evolution equations . Incorporating symmetries 

of DDs [|j imposes rather strong restrictions on realistic models of NFPDs. A possible strategy for a self-consistent 
model building is to use nonperturbative or phenomenological approaches (MIT bag, quark models, QCD sum rule 
ideas, etc.) to construct double distributions at low normalization point and then evolve them to higher values. 

The evolution equation for the nonsinglet quark double distribution was derived in Ref. where its analytic 
solution was also given. Evolution of the gluon distribution in pure gluodynamics was discussed in Refs. In 
this paper, we present also a full set of evolution equations for the flavor-singlet case and derive a solution following 
the method of Refs. An independent study of singlet evolution based on our approach was performed in 

Ref. JlO[ . Evolution equations for various versions of nonforward distributions can be found in 
convenient way to obtain the relevant evolution kernels is to use the universal light-ray evolution kernels 
The evolution of nonforward distributions was studied numerically in refs. pT]JT^ ,P, p^ , p^ . 

In the present paper, we incorporate the spectral and symmetry properties of double distributions to construct 
some simple models for DDs. Using the relations between DDs and NFPDs/OFPDs, we derive models for the 
latter and show that using the formalism of double distributions we can easily explain characteristic qualitative 
and quantitative features of the evolution of nonforward distributions observed in Refs. pl] , p7t . 




II. BASIC DEFINITIONS 



The kinematics of the amplitudes of the DVCS process "y*{q)N{p) j{q')N{p') and hard electroproduction 
^*{q)N{p) M{q')N{p') can be specified by the initial nucleon momentum p, the momentum transfer r = p — p' 
and the momentum q' of the final photon or meson. To get a Bjorken-type scaling limit, one should also keep the 
invariant momentum transfer t = small compared to the virtuality —Q^ = (q — r)^ of the initial photon and 
the energy invariant p ■ q = nipV. The essential features of the hard electroproduction processes (DVCS included) 
can be most easily demonstrated if we set q'^ = 0, = 0, = and use p, q' as the basic light-cone (Sudakov) 
4- vectors. It is easy to see that the requirement p'^ = [p + r)^ = p^ reduces in this limit to the condition p ■ r — Q 
which can be satisfied only if the two lightlike momenta p and r are proportional to each other: r = C,p, where C 
coincides with the Bjorken variable C, = xsj = /"^{p ■ q)- The latter satisfies the constraint < xsj < 1. For 
small but finite t and TOp, the momentum transfer r still must have a non-zero plus component r"*" = Cp^- It also 
may have a transverse component rj_. 

In the pQCD factorization treatment of hard electroproduction processes, the nonperturbative information is 
accumulated in the nonforward matrix element {p ~ r\ ip{0)ip{z) \ p) (we use here as a generic notation for quark 
{ijj) or gluonic (G) fields). It depends on the relative coordinate z through three invariant variables (pz), (rz) and 
z^. In the forward case, when r = 0, one gets the usual parton distributions by Fourier transforming the light-cone 
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projected (i.e., = 0) matrix element with respect to (j>z). In the nonforward case, we can try to start with the 
general Fourier representation 

/• OO POO /• oo 

(p-r|(^(0)(/3(z)|j5) =X((pz),(rz),z2;i,m2) = / dx\ dy e-"^P'''>''y^'''-'> p{x,y,iy;t) e-'''''di^ (2.1) 



oo — oo — oo 



with respect to all three z-dependent invariants. The Fourier transform p[x,y^v\t) can be called a triple distri- 
bution. Note that the generous (—00,00) limits for all three variables x,y,v serve for a most general function of 
(zp), (zr) and z^. However, incorporating information that the Fourier transformation is written for a function M 
given by Feynman integrals having specific causality properties, one arrives at more narrow limits: v runs from 
to 00, a: is between —1 and 1 while y is between and 1 (this was proven in |^ for any Feynman diagram using the 
approach of Ref. [^). To interpret the x- variable as the fraction of the initial momentum p carried by the relevant 
parton, it makes sense to separate integration over positive and negative x components and redefine x —x and 
y 1 — y for the negative x component. After that, the a;- variable is always positive and x and y are further 
constrained by inequality < x + y < 1 These spectral conditions can be summarized by the following 

representation 

{p~r\^{OMz)\p)= (2.2) 

fi i-i r 



JQ Jo 



-'^(P-)-'y(^-)^x, y, V- 1) + e"(P^)-*(i-«)('-^)$(x, y, v- 1) 



9{x + y < 1) dx dy , 



in which <i>(x, y, v. t) and $(2;, y, v\ t) result from positive- a: and negative- a; components of p{x, y, v] t) respectively. 
In particular, for quark operators, ^(a;, y, v; t) can be interpreted as the quark distribution while ^{x, y, v; t) as the 
antiquark one (a more detailed discussion is given in the next section). Similarly, y and [1 — y) can be interpreted as 
the fractions in which the momentum transfer r is shared among the two fields of the composite operator ip[Q)(p[z). 
Finally, the v variable characterizes the virtuality of these fields. For a light-cone dominated process, the leading 
term is given by the z^ — > limit of the nonforward matrix element, i.e. by zeroth moment of '^{x,y,v]t) with 
respect to v 



/>oo 

F{x,y;t) = / $(x, 

JO 



y,v\t)dv, (2.3) 



where F{x, y; t) is the double distribution. 

For a lightlike interval = 0, one can treat z as having only light-cone "minus" component, and then the scalar 
products {pz), (rz) project out the "plus" components of general (non-lightlike) momenta p and r. This allows to 
give a parton interpretation of F(x,y]t) as a probability amplitude for the active parton to carry fractions xp^ 
and yr~^ of the plus components of the external momenta r and p. Though the momenta p^ and r+ can be treated 
as proportional to each other r^ = Cp^j and r"*" specify the "-I-" -momentum flow in two different channels. 
For r+ = 0, the net "-I-" -momentum flows only in the s-channel and the total "-I-" -momentum entering into the 
composite operator vertex is zero. In this case, the matrix element is analogous to a distribution function. The 
partons entering the composite vertex then carry the fractions Xip'^ of the initial proton momentum (—1 < Xi < 1). 
When Xi is negative, we interpret the parton as belonging to the final state to secure that the integral always runs 
over the segment < x < 1. In this parton picture, the spectators take the remaining momentum (1 — x)p'^ . On 
the other hand, if the total "-I-" -momentum flowing through the composite vertex is r+, the matrix element has the 
structure of a distribution amplitude in which the momentum r+ splits into the fractions yr+ and (1 — y)r+ = yr+ 
carried by the two fields that appear in the vertex. In a combined situation, when both p^ and r+ are nonzero, the 
initial parton takes xp'^ + y"*", while the final one carries the momentum xp'^ — yr^ . For r = 0, we get the forward 
matrix element which is parametrized by the usual parton dist ribu t ions f{x). This gives reduction relations [0,00] 



connecting double distributions with the usual ones (see Eqs.(3.4),( |3.5| ) below) 



III. QUARK AND GLUON DISTRIBUTIONS 

For quark operators, the double distributions are defined by the following representation 

(p', s' I M'^)zE{^. A)Mz) \P.s)\ ,2=0 (3.1) 
= u{p\s')zu{p,s) ^ (e-^-(P-)-*^('-^)i^,(a;,y;t)-e"(P^)-*^('-^)Fa(a;,y;t)) 0{x + y<l)dxdy 

+ ^ u{p\ s'){zf - fz)u{p, s) J^' (e-"(P^)-^«(''^)Xa(x, y; t) - e'^^^'^-'y^'-'^ K-^{x, y; t)) e{x + y < 1) dx dy 
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for parton helicity-averaged ones and by 



ip', s' I M0)h5Ei0, z; A)M^) \P,s)\ ,2=0 (3.2) 



1 /.I 



u{p', s')z^5u{p, s) 1^ (e-"(f-)-*«('-^)G„(a;, y- 1) + e-(f^)-*^(-)Ga(a:, y; t)) e{x + y< 1) dx dy 



+ ^uip',s'h5u{P,s) 11 (e-"(^'^)-^^('-^)Pa(x,j/;t)+e"(P-)-^«('--)Pa(a;,y;t)) eix + y<l)dxdy 



1 <.i 



JQ 

in the parton helicity-sensitive case. Here and in what follows we adhere to the "bar" convention y — l — y,x— 1 — x, 
etc., for momentum fractions and use the "Russian hat" notation ^az" = z. As usual, u{p' , s'), u{p, s) are the Dirac 
spinors for the nucleon. In this definition, we explicitly separate quark and antiquark components of the double 
distribution. Note that such a separation is unambiguous: in the Fourier representation, it is completely determined 
by the sign of the x-term in the exponential. 

To clarify the physical meaning of separating the DDs into two components, it is instructive to consider the 
forward limit r = in which the matrix element is parametrized by usual parton densities, e.g., in the helicity 
averaged case 

{p I Via(0)zi?(0, z; A)M^) I p) I ,.=0 = u{p)zu{p) J^' (e-"(f^)/a(a:) - e"(f^)/a(x)) dx . (3.3) 

The exponential factors accompanying the quark and antiquark distributions reflect the fact that the field ip{z) 
appearing in the operator ipiO) . . . i^iz) consists of the quark annihilation operator (quark with momentum xp 
comes into this point) and the antiquark creation operator (i.e., antiquark with momentum xp goes out of this 
point). To get the relative signs with which quark and antiquark distributions appear in these definitions, we should 
take into account that antiquark creation and annihilation operators appear in V'(O) • . . ip{z) in the opposite order. 
Comparing the expression ( |3.3| ) with the r — limit of the definitions for DDs, we obtain "reduction formulas" 
relating the two components of the double distributions to the quark and antiquark parton densities, respectively: 

Fa{x,y;t^O)dy = fa{x) ; / Fa{x,y;t ^ 0) dy ^ fa{x) , (3.4) 

/o Jo 

and similarly for the helicity-sensitive case: 

l—x pl—x 

Ga{x,y;t = 0)dy = Afa{x) ; / Ga{x,y;t = 0) dy ^ Afa{x) . (3.5) 

JQ 

The reduction formulas tell us that integrating the double distribution Fa{x, y;t — 0) over a vertical line x =const in 
the {x, y)-plane, one gets the quark density /a (a;) while integrating its counterpart Fa{x, y;t = 0) gives the antiquark 
density fa{x). This is an illustration of our statement that Fa{x,y) and Fa{x,y) are independent functions. In 
particular, Fa{x,y ) co nt ains the valence component (reducing to /™'(a;)) absent in Fa{x,y). 



Our definitions (3T), ( |3.2[ ) reflect the results of the a- representation analysis ||5| that the plus component of the 
momentum of the particle (either quark or antiquark) going out of the hadronic blob can be written as xp^ + yr^ 
with both X and y positive and x + y < 1. This is in full compliance with the parton model based expectation 
that the initial hadron splits into an active parton and spectators which both carry positive fractions of its plus 
momentum. To show the positivity of the plus momentum component for spectators, we should explicitly take 
into account that, in the kinematics of DVCS and hard electroproduction processes, the plus component of the 
momentum transfer r = p ~ p' is positive r+ — C,p'^ > 0. Requiring that the plus component of the final hadron 
momentum is also positive, we conclude that < C < 1- Hence, < x + yC < 1 (since < C < 1 and < x + y < 1), 
i.e., the plus component of the momentum carried by spectators is also positive. On the other hand, the parton 
"going back" has the momentum whose plus component xp^ — yr^ = {x — yC)p^ may be either positive or negative, 
depending on the relationship between x, y and C- When (x — yQ is negative, one may wish to interpret such a 
parton as an antiparton leaving the hadron together with the initial parton. One should remember, however, that 
the double distributions F{x, y; t) "know nothing" about the magnitude of the skewedness C,: they are universal 
functions describing flux of p^ and r+ independently of what the ratio jp'^ might be. As we explained above, 
the quark DDs are unambiguously divided into two separate components Fa{x, y\ t) and Fa{x, y; t), but there is no 
further subdivision inside them based on interrelation between the values of x and y. 

In a similar way, we can introduce double distributions for the gluons 
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u{p\ s')zu{p, s) (z ■ p) 



1 /.I 



^0 



-ix{pz)—iy(rz) i ix(pz) — iy{rz) 



(3.6) 

{x + y < 1) Fg{x, y; t) xdxdy + "Kg" — term. 



{p',s' I z^z,Gl^iO)Eab{0,z;A)Gi,{z) \ p,s)\,2=o 



u{p\ s')zj5u{p, s) {z ■ p) 



-ix{pz)—iy{rz) ix{pz)—iy{rz) 



(3.7) 

{x + y < 1) Gg{x,y;t) xdx dy + "Pg"— term. 



There are no "antigluons" , so the positive-x and negative- a; parts are described by the same function. Note that 
our definition of the gluon double distributions here difi^ers from that used in our earher papers [^|J^,p^ by an extra 
factor of X in its right hand side. This form is more convenient for appUcations of the method of Refs. |^,^,^ to 
solve evolution equations for double distributions in the singlet case. The choice made above corresponds also to 
the simplest form of the reduction formulas 



Fg{x,y;t = 0)dy ^ fg{x) 



Ggix,y;t^O)dy = Afgix). 



(3.8) 



Another amb iguity in the definition of the gluon double distribution is related to the overall factor {z ■ p) in the 
rhs of Eqs.(^.6|) and ( ^.7[ ). Instead of it, we could take, e.g., (z • p') or (z • r) (such a choice is utterly inconvenient 
for taking the forward limit) or (z • P) where P = {p + p')/2 is a symmetric combination of the initial and final 
momenta. The latter choice (made in r ef. |9|| ) is more convenient for the studies of symmetry properties of the 
gluon DDs. Our choice made in Eqs. (p. 6]) and ( p.7[ ) simplifies the expressions for off-diagonal (QG and GQ) 
evolution kernels (see Eq.( 4.12| ) below). 

The flavor-singlet quark operators 



Oq{uz,vz) = X! 9 



'ipaiuz)zE{uz, vz; A)iljai'vz) — ipaivz)zE{vz, uz; A)ipa(uz) 



(3.9) 



and 



AOq{uz,vz)^J2^ 



■ipa{uz)z"/5E{uz, vz; A)ipa{vz) + ■>Jja{vz)z-f5E{vz, uz; A)')pa{uz) 



(3.10) 



are expressed in terms of double distributions Fq{x, y; t), Gq{x, y; t), etc. specified by 
{p\s'\Oq{uz,vz) I 1.2^0 = u{p',s')zu{p,s) (e— 

_^^xvipz)-^yvirz)-^^u{pz)-^u{rz)^ Fq{x , y; t) e{x + y < 1) dx dy + "Kg" -term, (3.11) 



{p',s' \ AOq{uz,vz) \ p,s)\ ^2^0 = u{p',s')z-f5u{p,s) I J ^ (e' 



-ixv{pz) — iyv(rz)-\-ixu{pz) — iyu{rz) 



^^^xv{pz)-^yv{rz)-^xuipz)-^u(rz)^^J Q q{x , y; t) 9{x + y < 1) dx dy + "Pq" -teil 



(3.12) 



They are given by the sum of "a -|- a" distributions: 

Nf 

FQ{x,y;t) ^Y,^Fa{x,y;t) + Fa{x,y;t)) 

a=l 



N, 



Gq{x, y; t) = Y.^Ga{x, y; t) + Ga{x, y; t)) . (3.13) 
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IV. EVOLUTION EQUATIONS 



The QCD perturbative expansion for the matrix element in Eq.(2T) generates Inz^ terms. As a result, limit 
^ is singular and the distributions F{x,y]i), etc., contain logarithmic ultraviolet divergences which require 

an additional i?-operation characterized by some subtraction scale /Lt: F{x^ y; t) F{x, y; 1 1 /x). The /z-dependence 

of F{x^ y;t \ ^) is governed by the evolution equation 



i^F,{x,y;t\fi)^ ^ Y^R-''{x,y;^,rj)Fi 



(4.1) 



where a,b — G,Q. A similar set of equations, with kernels denoted by AR°'''{x,y;£^,ri) prescribes the evolution of 
the parton helicity sensitive distributions G'°(x, y;t\ fi). Since the evolution kernels do not depend on t, from now 
on we will drop the i-variable from the arguments of F{x, y; 1 1 /i) in all cases when this dependence is inessential 
(likewise, the /i- variable will be ignored in our notation when it is not important). 

Since integration over y converts Fa (x, y; t — \ fi) into the parton distribution function fa (x \ /i) , whose evolution 
is described by the DGLAP equations [|§-|D 



M|^/a(x|M) 



(4.2) 



the kernels R {x,y;^,ri; g) must satisfy the reduction relation 







R-''{x,y;^,Tj;g)dy=^P-\x/C,g). 



(4.3) 



Aternatively, integration over x converts Fa{x, y]t — Q \ ^) into an object similar to a meson distribution amplitude 
(DA), so one may expect that the result of integration of i?"''(x, y; ^, 7y; g) over x should be related to the kernels 
governing the DA evolution |23J24]. For the diagonal kernels the relations are rather simple: 



R'^'^{x,y-^,ri-g)dx^V'^'i{y,r^-g) 



(4.4) 



for the quark kernel and a slightly more complicated expression for the gluon kernel: 



i?^^(x, y- ^, ?7; g)dx = (y, r/; g) 



(4.5) 



The x/^ factor appears because of the extra x which was added in the definition of the gluon DD by analogy 
with the definition for the usual gluon densities. The nondiagonal kernels R'-^^ and R^'~^ obey more complicated 
reduction formulas (see the Appendix). 

The reduction properties of the diagonal evolution kernels can be illustrated using the explicit form of the 
QQ-kernel: 



+ 



x,y;^,Vi9) 
0{O < xli < l)x/^ 



9(0 <x/^< mm{y/r], y/r]}) - -6(1 ~ x/^)S{y - -q) 



(4.6) 



-5{x/i - y/i]) + ^S{x/^ - y/rj) 



26{i-x/05{y-v) 



■ dz 



Here the last (formally divergent) term, as usual, provides the regularization for the l/(a; — ^) singularities present in 
the kernel. This singularity can be also written as l/{r] — y) for the term containing 6{x/^ — y/r]) and as 1/(77 — y) for 
the term with d{x/^ — y/fj). Depending on the chosen form of the singularity, incorporating the 1/(1 — z) term into a 
plus-type distribution, one should treat z as a;/^, y/rj or y/fj. One can check that integrating R^^{x, y',£,,Ti; g) over 
y 01 X gives the DGLAP splitting function P'^'^ {x / g) and the DA evolution kernel {y^rj\ g) , respectively: 



pQQ[z,g)^^cJ\ 



1 - z 



V<^<^[y,,;g)^^C,[(y- 



V - y 



0{y < r,) 



Oiy > v) 



(4.7) 
(4.8) 
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Here, "+" denotes the standard "plus" regularization 

A convenient way to get explicit expressions for ^, 77; is to extract them from the kernels w) 

describing the evolution equations for the light-ray operators p4[p^p^,pT 



H-^Oa{Q,z)^ [ [ B'''>(u,v)Ob(uz,vz)0(u + v<l)dudv. (4.9) 
Jo Jo V 

Since the definitions of the gluon distributions F^(x,y;t), G^{x,y]t) contain an extra (pz) factor on the right- 
hand side, which results in the differentiation d/dx of the relevant kernel, it is convenient to proceed in two steps. 
First, we introduce the auxiliary kernels r°'''{x,y;£,,ri;g) directly related by 

r°''(x,?/;^,77;g) = / [ 5{x -£_{!- u - v)) S{y - u - - u - v)) B''^{u,v) 9{u + v < 1) dudv 
Jo Jo 

= lB'^\y-rjx/^,y-nx/0 (4.10) 
to the light-ray evolution kernels B°'''{u,v). The second step is to get the _R- kernels using the relations 



RQ(^{x,y;C,ri;9) =r^^{x,y;^,rj;g) , R^^{x,y;i,mg) = -r^^{x,y;^,r];g) (4.11) 

X 



d_ 

dx 



(^xR'^^{x,y;tv;9?j ^ -r"^^ {x,y; ^,r]; g) , R^"^ {x,y; ^,7j; g) = ^r^"^ {x,y;^,7^; g) . 



(4.12) 



Hence, to obtain R'^'^{x,y;£_,ri;g), we should integrate r'-^'^{x,y;^,r];g) with respect to x. We fix the integration 
ambiguity by the requirement that R'^'^{x, y; ^, 77; g) vanishes for x > 1. Then 



1 



R'"^{x,y;^,V;g) = - / r'^'^ {x,y; 9) di . (4.13) 



This convention guarantees a simple relation (4.3) to the DGLAP kernels. Explicit expressions for the evolution 
kernels and discussion of evolution equations in the singlet case is given in the Appendix (see also ref. pO[). 



V. PARTON INTERPRETATION AND MODELS FOR DOUBLE DISTRIBUTIONS 



The structure of the integrals relating double distributions with the usual ones 

fa.a,g{x)^ / Fa,a,g{x , y)dy (5.1) 

Jo 

[where F{x,y) = F(x,y;t — 0) ] has a simple graphical illustration (see Fig.||a). The DDs F{x,y) live on the 
triangle defined by < a;, ?/, a; -I- y < 1. Integrating F{x,y) over a line parallel to the y-axis, we get f{x). The 
reduction formulas and the interpretation of the x-variable of F{x,y) as a fraction of the momentum suggests 
that the profile of F{x,y) in the x-direction is basically driven by the shape of f{x). On the other hand, the 
profile in the j/-direction characterizes the spread of momentum induced by the momentum transfer r"*". Hence, 
the y-dependence of F{x,y) for fixed x should be similar to that of a distribution amplitude ipijj). By analogy 
with, e.g., the pion distribution amplitude ifj^^y), which is symmetric with respect to the change y 1 — y, one 
may expect that the distribution of the r-momentum between the two partons described by the same field should 
also have some symmetry. However, the symmetry cannot be as simple as ?/ <-> 1 — y since the initial p and the 
final p' = p — r momenta arc not treated symmetrically in our description: the variable x specifies the fraction of 
the initial momentum p both for the outgoing {xp + yr) and incoming [xp — (1 — y)r) partons. To treat p and 
p' symmetrically, we should interpret x for the returning parton as the fraction of the final hadron momentum 
p' = p — r, i.e., rewrite its momentum xp — [1 — y)r as x{p — r) — [1 — x ~ y)r. Hence, the symmetry of a double 
distribution F{x^ y) may be only with respect to the interchange ?;<->l — x — 

Another way to make the symmetry between the initial and final hadrons more explicit is to use P = {p + p')/2 
and r as the basic momenta rather than p and r (cf. |0,^|J^ ) writing the momenta of the partons as xP + yr and 
xP — (1 — y)r. Then the y 1 — x — y symmetry corresponds to y <-> 1 — y symmetry. The variable y changes in 
the interval x/2 < y < (I — x/2). Writing y as y ~ {I + a) /2, we introduce a new variable a satisfying a symmetric 



7 



constraint —x < a < x, where x = 1 — x. The y <-> I — x — y symmetry now converts into a <-> —a symmetry. 
Finahy, rescahng a as a = xf3 produces the variable (3 with a;-independent hmits: —1 < /3 < 1. Written in terms 
of x and /3, a modified double distribution F(x, (3) obeys the reduction formula 



1 

F{x,(3)dp = f{x). (5.2) 



It is instructive to study some simple models allowing to satisfy this relation. Namely, let us assume that the profile 
in /3-direction is a universal function g{[3) for all x, i.e., take the factorized ansatz 

F{x,P)^-^f{x)g(P), (5.3) 
1 — a; 

with g{[3) normalized by 

fi 

5(/?)d/3 = l. (5.4) 

Possible simple choices for g{j3) may be 5{(3) (no spread in /3-direction) , |(1— /3^) (characteristic shape for asymptotic 
limit of quark distribution amplitudes), j|(l — /3^)^ (asymptotic shape of gluon distribution amplitudes), etc. In 
our original variables x,y, the factorized ansatz can be written as 

F{x,y)=^^^f{x), (5.5) 

where h(x,y) is a function symmetric with respect to the interchange j/<->l — a; — y. A trivial observation is that 
the variable x itself is given by a combination [1 — {1 — x — y) ~ y] symmetric with respect to the y <-» 1 ^ x — y 
transformation. The normalization function h{x) is specified by 

h{x) = / h{x,y)dy. (5.6) 



For the three simple choices mentioned above, the model (5^) gives 

F('^Hx,y) = 6{y~x/2)fix) , F(^\x,y) ^ ^^^(^^I^^^^^ ' ^^'^^^^v) = ^^^'[IZ"^)^ /(^) ' (^.7) 
In a similar way, one can construct ansatze for functions F{x,y]t) involving nonzero t values. 



VI. RELATION TO NONFORWARD DISTRIBUTIONS 



The nonforward matrix elements accumulate process-independent information and, hence, have a quite general 
nature. The coefficient of proportionality between p"*" and characterizes the skewedn ess of mat rix elements. The 
characteristic feature implied by representations for double distributions (see, e.g., Eqs.( |3.l|) , (^^)) is the absence of 
the (^-dependence in the DDs F{x, y) and G(x, y). An alternative way to parametrize nonforward matrix elements 
of light-cone operators is to use the ratio C, = jp^ and the total momentum fraction X = x -\-yC, independent 
variables. Taking into account that for a lightcone dominated process only one direction for z gives the leading 
contribution, one can do the change (rz) = C{pz) directly in our definitions of double distributions. As a result, the 
variable y would appear there only in the x + y(^ = X combination, where X can be treated as the total fraction 
of the initial hadron momentum p carried by the active quark. If we require that the light-cone plus component of 
the final hadron momentum (i.e., p"*" — r"*") is positive, then < C < 1- Using the spectral property < x + y < lof 
double distributions we obtain that the variable X satisfies a similar "parton" constraint < A < 1. Integrating 
each particular double distribution Fa.a.g{X — y(, y) over y gives the nonforward parton distributions 

J^^'^'a{X)^9{X>0 FaraAX~yC,y)dy + e{X <C) Fa.-a.g{X - yCy) dy , (6.1) 

Jo Jo 

where ^ = 1 ~ (. The two components of NFPDs correspond to positive {X > () and negative {X < C,) values of 
the fraction X' = X — C, associated with the returning parton. As explained in refs. [HJ^], the second component can 
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be interpreted as the probability amplitude for the initial hadron with momentum p to split into the final hadron 
with momentum (1 — C,)p and the two-parton state with total momentum r = C,p shared by the partons in fractions 
Yr and (1 - Y)r, where Y = X/C, (see Fig.|). 

For the gluon DDs, the y <-> \ — x — y symmetry holds only if, instead of {z ■ p), one uses the symmetric 
overall factor (z • P) in the definitions (3.6),(3.7)[]. The use of such a definition of the gluon DDs is implied in 
this section. Furthermore, the nonforward gluon distribution J-^{X) is obtained by integrating xFg{x,y)\x=x-yC- 



To simplify notations, it will be also implied below that, for the gluons, F{X — yC,y) in Eq.( |6.l[ ) corresponds to 
{X-yC)Fg{X-yC,y). 

The basic distinction between double distributions F{x,y) and nonforward distributions JF^(X) is that NFPDs 
explicitly depend on the skewedness parameter They form families of functions J-'^'"''^ (X) whose shape changes 
when C is changed. The fact that the functions ^({X) corresponding to different C's are obtained by integrating 
the same double distribution F{x,y) imposes essential restrictions on possible shapes of J-'({X) and on how they 
change with changing C. The relation between NFPDs and DDs has a simple graphical illustration on the "DD-life" 
triangle defined hy < x,y, x + y < 1 (see Fig.|| ). To get T(^{X), one should integrate F{x,y) over y along a 
straight line specified hy x = X — C,y. Fixing some value of C, one deals with a set of parallel lines corresponding to 
different values of X. Evidently, each such line intersects the x-axis at x = X . The upper limit of the y-integration 
is determined by intersection of this line either with the line x + y = 1 (this happens ii X > Q or with the y-axis 
(if X < C). The line corresponding to X = C separates the triangle into two parts generating two components of 
the nonforward parton distribution. In the forward case, when C = 0) there is only one component, and the usual 
parton densities f{x) are produced by integrating F{x,y) along the vertical lines x =const (see Fig.^j). In case 
when X > C,, looking at the integration line for the nonforward parton distribution T(;{X) one can see (Fig.||b) that 
it is inside the space between the integration lines giving the usual parton densities f{X) and f{X') corresponding 
to the momentum fractions X, X' = X — ^ of the initial and final parton. Assuming a monotonic decrease of 
the double distribution F{x, y) in the a;-direction and a universal profile in the y-direction, one may expect that 
Tc^{X) is larger than f{X) but smaller than f{X'). Inequalities between forward and nonforward distributions 
were recently discussed in refs. |l7| , ^5|j26| . They are based on the application of the Cauchy- Schwartz inequality 



s 



H{py,x'p,s\H{py,xp,s)\^ <Y,{Hip),xp,s\Hipy,xp,s)J2{Hip');x'p,s'\Hipy,x'p,s'), (6.2) 



to the nonforward distributions J^^(X) written generically as 

s 

where \ H{p); Xp, S) describes the probability amplitude that the hadron with momentum p converts into a parton 
with momentum Xp and spectators S. The forward matrix elements arc identified with the usual parton densities 

Y,{H{p);Xp,S\H{p);Xp,S) -/(X). (6.3) 

s 

Notice that the hadron momentum in the second forward matrix element is p' = Qp, hence the argument of the 
relevant parton density is X'/C, and one has 

Y,{H{py,X'p,S\H{py,X'p,S) = Y,{H{py,X'p' /lS\H{py,X'p' lis) ^ f{X' lOIC. (6.4) 

s s 

As a result, we obtain (compare with p5|j2^) 



rl{x)<^f{x)f{x'/oic< 



2yw 



/(x) + /(x7c) 



(6.5) 



In other words, the functions involved in the bound for !Fq{X) are /(X) and /(X2) where the fraction X2 = X'/C 
is larger than X' |23,E0]. One can see that X2 is given exactly by the a;- value of the intersection point in which 



I am grateful to G. Filler for attracting my attention to this point. 
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the integration hne x = X — C,y giving the nonforward distribution !Fq{X) crosses the boundary hne a; + y = 1 
(see Fig.^). For the gluon nonforward distributions, one should take into account extra factors {zp), {zp') present 
in the definitions of forward distributions and the overall factor chosen in the rhs of definitions of the nonforward 
gluon distributions (see Eqs. (3J), ( |3.7|) ). If one uses the p ^ p' symmetric combination (z • P), then 



■g (symm) 



c 



{X)< 



1 



l-C/2 



f{x)f{x'/o< 



1 



2(1 -C/2) 



.f{x) + f{x'/o 



(6.6) 



It is clear that the whole construction makes sense only \i X' > Q (or X > Q. If X' < 0, the nonforward 
distribution corresponds to matrix elements {H{p'); Xp, X'p, S \ H{p), S ) which have no obvious relation to the 
usual parton densities. Furthermore, in our graph of Figjja, the left end of the line x = X — C,y in this case 
corresponds to x = 0, where the usual parton densities are infinite, and the inequalities become trivial. In fact, 
they are trivial even for the border point X = Another deficiency of the Cauchy-Schwartz-type inequalities is 
that they do not give the lower bound for nonforward distributions though our graphical interpretation suggests 
that T^iX) for AT > C is larger than /(AT) if the x-dependence of the double distribution F{x,y) along the lines 
y — kx is monotonic. 

To develop intuition about possible shapes of nonforward distributions, it is instructive to derive the NFPDs corre- 
sponding to three simple models specified in the previous section. In particular, for the F^^\x, y) = 6{y — x/2)f(x) 
ansatz we get 



:Ff\x) 



9{X>C/2) JX-(:/2 



l-C/2 



/ 



l-C/2 



(6.7) 



i.e., NFPDs for non- zero C are obtained from the forw a rd d istribution f{X) = J^q{X) by a shift and rescaling. 
Note that the model ( |6.7| ) satisfies the inequalities (6.5), (6.6) in the region X > for any function f{x \ Qo) of the 
f{x I Qo) = Ax^"-{1 — x) type provided that a > and 6 > 0. Using the relations 



Hix,C,t) = il-C/2)TciX;t) ; .i 



^-C/2 
l-C/2 



c 



2-C 



(6.8) 



between our nonforward distributions and Ji's off-forward parton distributions (OFPDs) H{x,£^;t) ||l|;^, one can 
see that the delta-function ansatz gives the simplest ^-independent model 

H<^^\x,^;t^O)^f{x) 

for OFPDs at t = |. It is worth noting that the MIT bag model calculation |2^ did produce a set of OFPDs 
which are almost independent of An evident interpretation is that the model constructed in ref. ||27[| strongly 
suppresses the redistribution of the momentum transfer among the constituents which results in a very narrow 
spread of F{x,y) in the y-direction. Even if such a picture is physically correct for a low normalization point 
Qo ~ 500 MeV, evolution to higher values Q ^ 1 GeV widens the y-profile of F{x,y) and evolved OFPDs would 
change their shape with ^, as was explicitly demonstrated through a numerical calculation by Belitsky et al. p^ ]. 

The evolution of nonforward distributions Tq{X\Q) was recently studied in refs. |ll|Jl7| , p!q| . As a starting 
condition, the authors assume that, at some low scale Qo, the nonforward distributions Tq{X \ Q) for all C have the 
same universal shape coinciding with that of the usual (forward) densities f{X, Qo). This assumption corresponds 
to the ansatz F{x, y \ Qq) = 5{y)f{x \ Qq) with double distribution being nonzero on the a;-axis only. This ansatz is 
not realistic, since it has no symmetry with respect to the y (l — x — y) interchange. However, evolution equations 
are applicable to any distribution and, just due to its asymmetric profile, this unrealistic double distribution has a 
very distinctive evolution pattern reflecting the restoration of the y l — x — y symmetry. Namely, the asymptotic 
functions F{x,y \ Q — > oo) are y {l — x — y) symmetric. In particular, both in pure gluodynamics and in QCD, we 
have Fg{x,y\Q ^ oo) ~ y^(l — a; — y)^ (see Ref. Q and the Appendix). Hence, one may expect that the evolition 
of Fg[x, y I Q) shifts its crest towards the y = x/2 line and also makes the y-shape of the double distribution wider. 
To see whether the results of Refs. [ [ll|jl7| ] reflect this expectation, we introduce a general model with a narrow 
y-dependence: F^^\x,y) = 5{y — kx)f{x) (in what follows, it will be referred to as the "fc-delta ansatz"). This 
double distribution is concentrated on the y = kx line and gives 



■'■Since hadrons are massive, t = is outside the physical region; hence, the t limit should be understood in the sense 
of analytic continuation. 
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for nonforward distributions. In case of two other models, simple analytic results can be obtained only if we 
specify a model for f{x). For the "valence quark" -oriented ansatz F'^^">(x,y), the following choice of a normalized 
distribution 

is (a) close to phenomenological valence quark distributi ons a nd (b) produces a simple expression for the double 
distribution since the denominator (1 — x)^ factor in Eq.(6.10) is canceled. As a result, the integral in Eq.(6.1) is 
easily performed and we get 



T^'\x)^^lx'-^{CaX-2{X-C))+9{X>C)(^^^] {CaX + 2X0\. (6.11) 



2-a 



Resulting curves for J^^^^(X) with a = 0.5 and C = 0.05,0.1,0.2,0.4 are shown in Fig.^ A characteristic feature 
of each curve is a maximum located close to the relevant border point X = ( and slightly shifted to the left from 
it. Note that both the functions J-'^^\x) and their derivatives {d/dX)T^^\x) are continuous at X ^ (^. The 
latter property is secured by the fact that F(^)(x,2/) vanishes at the upper corner x = 0,y — I. The (1 — x)^ 
denominator factor for the "gluon-oriented" ansatz F^'^\x, y) is canceled if one takes the model f{x) ^ x^°'{l — x)^ 
which, fortunately, is also consistent with the a; — > 1 behavior of the phenomenological gluon distributions. It is 
well known that the values of nonforward distributions JF^(X) taken at the border point X — ( determine 
imaginary parts of DVCS and hard electroproduction amplitudes. An interesting question is the relation between 
the usual distributions /(C) and the values .?x(C) of nonforward distributions at the border point. It is easy to 
calculate that for the k = 1/2 delta ansatz F'-^^x, y) this ratio is given by 

i?(°)(C) ^ Tf\C)lf{0 = l^^lf/'^f^^ ■ (6-12) 

It is larger than 1 for any monotonically descreasing function f{x), i.e., the nonforward distribution Tq{C) in this 
case is larger than /(C). In the small-C limit, R^'^^Q is completely determined by the small-x behavior of f{x), 
and the expression for R^^\C,) simplifies to 

i?^°HC) lc«i « ^ • (6.13) 

Hence, if f{x) has a purely powerlike behavior f{x) ^ for small x, then R''^\C, ^ 0) = 2°(1 + 0(C)), i-c, for 
small Cj the ratio of the nonforward distribution !Fq{C,) and the usual parton density /(C) is practically constant, 
deviating from the C = limiting value by 0(C) terms only. The limiting value in this case is 1.41 for a — 0.5 and 
1.23 (1.15) for a = 0.3 (a = 0.2). However, if f{x) is a sum of two different powerlike terms Ax~°' + Bx^^ or if it 
contains logarithms, e.g., f{x) ~ x~°' ln(l/a;) for small x, then the C-dependence is more pronounced. In the latter 
case 

and there is a visible deviation from the limiting C ^ value for all accessible C^ on the ln(l/C) scale, the 
C-dependence of the ratio i?^°^(C) cannot be neglected even for C ~ 10^'^. 

For a general A:-delta model F^\x, y) = 5{y — kx)f{x), the ratio J-q{C,) / f{C) for small C can be approximated by 
/(C(l — k))/ f{Q which again gives a C-independent constant (1 — fc)^" for a purely powerlike function f{x) ~ a;~" 
while the ln(l/x)-factor would modify the constant by [1 -I- ln(l — fc)/lnC]. 

If one uses the "valence quark" -oriented ansatz F'-^^x, y) with a simple powerlike behavior J{x) ^ x^" for small 
x, the ratio is given by 

i?(^)(C) . TfiO/m - ^,_cni-'a/2)il-a/3) ' ^^'^^^ 
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Just like in the previous example, the nonforward distribution ^^(C) is larger than /(C) for all positive a. For 
small C, the ratio tends to 1/(1 — a/2)(l — a/3), e.g., to 1.6 for a — 0.5 which is the usual choice for valence 
quark distributions (for comparison, taking a = 0.3 (a = 0.4) gives 1.3 (1.44) for i?'-^-'(C 0)). For small a, this 
result can be translated into i?(°)(C 0) « e""/^ « /(e^/^C)//(C), which coincides with the ratio R'^°\C, 0) 
for the modified narrow ansatz F^\x,y) — 5{y — kx)f{x) with k « 0.56. Hence, for Tq{C,) the widening of the 
y-distribution can be approximated by a narrow distribution shifted from y — x/2 upwards to the y kx line. 
Again, a logarithm ln(l/x) in f(x) at small a: would induce a visible ^-dependence for the R^^\C,) ratio even for 
very small C- 

Switching to the "gluon-oriented" ansatz F'^'^\x,y) with a purely power behavior f{x) ^ x^^ for small x, we 
obtain a similar expression 

^'"«'"^c''K'/^":'- (l-C)3(l-„/3)a-a/4)(l-a/5) ' 

which is close to (2.17)" for small a. To approximate this result by the delta ansatz F^^\x,y), one should take 
k w 0.54. The effective shift upward is smaller in this case because F^'^\x,y) is more narrow in the y-direction 
than F'^^\x,y). 

Choosing a, we should take into account that the nonforward gluon distribution J-^(X) reduces to Xfg{X) in 
the C ^ limit |^,||. Hence, /(C) in the above formulas should be understood as CfgiO- Now, if we make an 
old-fashioned assumption that Xfg{X) tends to a constant as X — > 0, then a — and i?'^^^(C) tends to 1 at small 
C, i.e., the nonforward distribution coincides in the small-C limit with its forward counterpart CfgiO- To 

get a more realistic gluon distribution Xfg{X) growing at small X one should use a positive parameter a. Taking 
a = 0.3, we get R'^^\( 0) « 1.27, and i?(2)(( ^ Q) « 1.17 (1.39) for a = 0.2 (a = 0.4)). 

These estimates for the ratio T^{X \ Q) / X fg{X \ Q) are close to those obtained in refs. | |ll| , [l7t where the nonfor- 
ward distributions {X \ Q) at high normalization point Q were constructed by applying evolution equations to an 
initial low normalization point Qo ansatz J-^iX \ Qq) which was assumed to have a universal C- independent shape 
coinciding with the usual distribution X fg{X \ Qo). In particular, Martin and Ryskin considered the evolution of 
the gluon NFPD in pure gluodynamics. They took Qq — 1.5 GeV'^ (two other choices Qq — 0.4 GeV^ and = 4 
GeV^ were also considered) and then evolved J^^{X \ Q) to higher values = 4,20, and 100 GeV^. They 

found that i?(10~^) « 1.3 for Q^ — 100 GeV^, which corresponds to a « 0.3 in our F^^^ model. This value is close 
to those used in phenomenologicalparametrizations of the gluon distributions. It should be also noted that the 
results for -R(C) obtained in Ref. [|7| have a nonnegligible C-dependence. This feature can be expected since the 
GRV gluon distribution which they use can be rather well approximated at Q^ = 4 GeV^ by a simple formula 

x/f^^(x,Q2 =4GeV2) « ia;-"-3ln(l/a;) 
y 4 

which works with 10% accuracy for x ranging from 10^^ to 10^^. In the pure gluodynamics approximation used 
in Ref. ||l^, its shape does not drastically change when evolved either to Q^ — 1.5 GeV^ or to Q^ — 20 and 100 
GeV2. 

As discussed above, the assumption that the nonforward distributions J-^ {X \ Qq) have a universal C-independent 

shape corresponds to the ansatz F^^\x, y \ Qo) — S(jj) f{x \ Qo), i.e., to the fc-delta ansatz with the vanishing slope 
k = 0. Modeling the evolved double distributions by a /c-delta ansatz with nonzero fc, we expect that, due to the 
restoration of the y 1 — x — y symmetry, the effective slope parameter k should increase with Q^. Namely, for the 
/c-delta ansatz, the ratio of the nonforward distribution {X) and the forward parton distribution f{x) = Xfg{X) 
is given by 

^c(-^) _ f{x-kcc/{i-kc)) 
^) = xf^ix) - (1 - kOfix) ■ ^"-''^ 

Taking f{x) = ia;-0-3ln(l/a;) and the Q^.^ependent slope k{Q^) = 0.3; 0.4; 0.48 for ^ 4-20 and 100 GeV^, 
respectively, we were able to reproduce the results of Ref. |l^] for a wide range of C parameters: C = 10~^, 10~^, lO""* 
and 10~^. The relevant curves, coinciding with those of Ref. within a few per cent accuracy, are shown in 
Fig.^. Hence, the increase of the ratio R{X, C) with observed in Refs. 1^,0 basically reflects the shift of the 
gluon double distribution from the x axis y = towards the symmetry line y = x/2. This effect, being an artifact 
of the initial conditions, plays the dominant role up to Q^ ~ 100 GeV^. As argued above, the C-dependence of the 
ratio may be traced to the fact that the gluon distribution xf^^^ (x) differs from a simple power x"". 



12 



Since the assumption Ti^{X \Qo) = f{X\Qo) is equivalent to the ansatz y | Qo) = ^{y)f{x\Qo) which is 
not symmetric with respect to the y ^ (1 — x — y) interchange, one should avoid using it as a starting condition 
for evolution. As explained earlier, a more realistic set of nonforward distributions 



e{x>c/2)_^fx-ci2 



l-C/2 



l-C/2 



(6.18) 



is generated by the F^^\x, y \ Qq) = 5{y — x/2)f{x) ansatz for the double distribution corresponding to skewedness- 
independent set of Ji's off-forward distributions. Comparing thes e two sets, one may be tempted to argue that for 
extremely small C, considered in Ref. jl^, C/2 terms in Eg. ( 6. 18 ) are inessential. Of course, C/2 can be neglected 
when subtracted from 1. However, for the X-values close to the border point X = C,, the shift by C/ 2 produces 
visible changes for functions having the X~°' behavior with a ~ 0.3. In the case of the ansatz ( |6.18 ), the ratio 



C I Q) = !Fq{X I Q) / X fg{X I Q) differs from 1 for all Q. For small C, the difference is significant only for X 
close to C- 

When a narrow double distribution has its crest on the y = x/2 line from the very start, there are no effects 
due to the shift of the crest, and the Q-evolution of R{X, C\Q) in the region X > C, reflects only the widening of 
the double distribution in the ?/-direction and the change of its profile in the x-direction. As we have seen, the 
widening of the double distribution changes the effective slope /c by a small amount only. Hence, for small ^ one 
can use the approximate formula 



THX I Q)|c«i « (X - CmfgiX - C/2 I Q) 



(6.19) 



for evolved distributions as well. In other words, the ratio R{X, ( \ Q) ioi X > ( and small C, can be estimated from 
existing results for the usu al glu on density f{X) = Xf c{X \ Q). 

Comparing the formula ( 6. IE ) with the relation (6.8) between our nonforward and Ji's off-forward distributions, 
one can conclude that Eq.( 6.19 ) is equivalent to a statement that at small ^ and x > ^ one can neglect the ^- 
dependence of the o ff-for ward distributions H{x; Again, such a statement is only nontrivial if a; ~ ^. To analyze 
the accuracy of eq.( |6.ig| ), we will construct an expansion of H{x; ^) in powers of ^. To this end, it is convenient to 
use the parton picture based on modified double distribution F(x, a) in which the plus component of the parton 
momenta is measured in units of that of the average hadron momentum P = {p + p')/2. The parton momenta then 
are xP+ (1 -|-a)7'/2 and xP— (1 — a)r/2 with a changing between —x and x. Defining r+/P+ = 2^ and x = x + ^a, 
one obtains the description in terms of the off-forward parton distributions H{x; ^) ||l|,^. The parton momenta are 
now {x + ^)P and {x ~ £,)P. In the region x > the OFPDs are obtained from F{x, a) by the integral 



(i-x)/(i-4) 



F{x — £,a, a) da . 



(6.20) 



Using the a —a symmetry of F{x, a), it is easy to see from this expression that the off-forward parton distribu- 
tions H(x; ^) are even functions of ^: 



3.21) 



Thi s resu lt was originally obtained by X. Ji [g5[ with the help of a different technique. Expanding the rhs of 
Eq.( |6.20 ) in powers of we get 



1 /-(i-*) d^F{i,a) 



2 



da 



dx 



(6.22) 



a — l — x- 
2 



where /(x) is the forward distribution. Hence, for small ^, the corrections are formally 0(^ ), i.e., they look very 
small. However, if F{x^a) has a singular behavior like a;~°, then 



d'^F{i,a) a{\ + a) 
dx"^ x"^ 



F{£,a) 



and the relative suppression of the first correction is 0{S^'^ /SP') i.e., the corrections are tiny for all x except for the 
region i ^ ^ where the correction has no parametric smallness. Nevertheless, even in this region it is suppressed 
numerically, because the a^ moment is rather small for a distribution concentrated in the small-a region. This 
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discussion shows that the formula ( 6.19| ) is not just an automatic consequence of the O(C^) nature of the first 
nonvanishing correction. It is easy to write expicitly all the terms which are not suppressed in the i ~ f — > limit 



^2k Ai-i) d'^''F{i,a) 



dx 



(6.23) 



The numerical suppression of higher terms is even stronger , an d the series converges rather fast. 
In terms of the off- forward distributions, the inequality (|6.5|) reads 



2yT^L + ? 



For the gluons, one should use the inequality (6.6), which leads to 



1 + C 



1-C 



(6.24) 



(6.25) 



Again, if one takes the model H{x,^) = f{x), the inequalities ( 6.24 ) and ( 3.25 ) are valid for any function f{x) of 
x'^il - x)^ type with a > 0, 6 > 0. 

So far we assumed in our models that DDs are finite everywhere on the "life triangle" . Consider, however, a 
situation when the partons emerge from a meson- like state (or glueball/pomeron in the gluon case) exchanged in 
the t channel. In this case, the partons just share the plus component of the momentum transfer r: information 
about the magnitude of the initial hadron momentum is lost if the exchanged particle can be described by a pole 
). Hence, the meson-exchange contribution to a double distribution is proportional to 5{x) 



propagator ~ 1/ (i — "m-M 
or its derivatives, e.g.: 



F^'\x,y) ^ 5{x) 



'M 



t 



(6.26) 



where tpM{y) is the distribution amplitude of the meson M. This contribution to the nonforward distribution is 
nonzero only in the < X < ( region: 



t) 



3.27) 



At the beginning, we described the nonforward matrix element of a quark operator by two functions F''(x,y) and 
F°'{x,y) corresponding to positive-x and negative- 2; parts of the general Fourier representation. Since x = for 
a meson-exchange contribution, it makes sense to treat it as a third independent component, i.e., to parametrize 
the nonforward matrix element by the sum ® ® F*^ . All three components contribute to the nonforward 
distributions in the < X < C region. However, the 8{x) terms do not contribute to the nonforward distributions 
in the > C region and to the usual parton densities f{x). For this reason, the 5{x) terms, if they exist, would 
lead to violation of sum rules (like energy-momentum sum rule) for the usual parton densities. 

Note that if the meson DA (/3(y) does not vanish at the end-points, the nonforward distribution does not vanish 
at AT = (the off-forward parton distributions H{x]^) in this case are discontinuous at x = ±^). As explained in 
ref. 1^, pQCD factorization for DVCS and other hard electroproduction processes fails in such a situation, because 
of the 1/A factors (1/(5; ± i^) factors if OFPD formalism is used) contained in hard amplitudes. It should be 
mentioned that a nearly discontinuous behavior of OFPDs for i = ± ^ was obtained in the chiral soliton model 
pof . Formally, the evolution to sufficiently high /i results in the functions vanishing at the end-point A = 0. A 
non-trivial question, however, is whether evolution starts at all in a situation when pQCD factorization fails. 



VII. SUMMARY 



In this paper, we duscussed the formalism of double distributions. We treated them as the starting objects in 
parametrization of nonforward matrix elements. An alternative description in terms of nonforward or off-forward 
parton distributions was obtained by an appropriate integration of the relevant DDs. Incorporating spectral and 
symmetry properties of double distributions, we proposed simple models producing self-consistent sets of non- 
forward distributions J^(^{X) and discussed their C-dependence and relation to usual (forward) parton densities. 



14 



Using a qualitative picture of the evolution of double distributions, we were able to explain and model the basic 
features of the evolution pattern of nonforward distributions observed in numerical evolution studies . In the 
Appendix, we present the set of evolution equations for double distributions in the singlet case and discuss their 
analytic solution. Work on numerical evolution of the nonforward distributions corresponding to realistic ansatze 
( |6.18 ) is in progress [Q. Another interesting problem for a future investigation is a numerical evolution of double 
distributions. 
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APPENDIX A: EVOLUTION EQUATIONS FOR THE SINGLET CASE 



As described in Sec. IV, the evolution kernels for double distributions can be conveniently obtained from the 
light-ray evolution kernels u). For the parton helicity averaged case, the latter were originally obtained in 

Refs. [p^lJTsl. Here we present them in the form given in Ref. H|: 



u, v) = ^Cf ( 1 + S{u)[v/v]+ + S{v)[u/u]+ - ^diu)div) ] , 



TT 

a. 



TT 
Os 

TT 
Os 

TT 



CF[2 + 5iu)5iv)j , 
Nf (1 + Auv - u-v) , 

Nc (4(1 + iuv-u-v) + 5{u)5{v) + { 5{u) 



Po 



2Nc 



Siv) 



dv 

V 



+ {u ^ v} 



(Al) 

(A2) 
(A3) 
(A4) 



As usual, /3o = 11 — |-/Vy is the lowest coefficient of the QCD /^-function. Evolution kernels for the parton helicity- 
sensitive case are given by | |l^ , [l3t 

(A5) 
(A6) 

(A7) 
(A8) 



(u,v) = B'^'^{u,v) 
AB^O (u, v)^^Cf (^6{u)6{v) - 2^ 

AB^'^iu, v)^—Nf{l-u--v), 



AB'"'"(u,u) = B''^(u,u) - 12 — Ncuv. 

TT 



At one loop, AR'^'^{x,y;£,,r];g) = R'^'^{x,y;£^,jy,g), and this kernel was already displayed in Eq. (4.6). Other 
kernels, including the i?*^*^(x, y; ^, ry; kernel originally obtained in Ref. are given by 



e{o<x/^<i)ix/o 
(1-^/0 



< x/^ < mm{y/ri, y/r]}) + S{1 - x/^)6{y - rj) 



f3o 
2Nc 



-6{x/£,- y/i]) + ^5{x/£,- y/rj) 
r] 7] 

1 



26{l-x/OS{y~v) 



dv 
l-v 



R^'-ix, y- C, ry; g) = AR^^{x, y; g) + l2 — N,-{y- rjx/Oiv - fjx/0 0(0 < x/^ < mm{y/7], y/fj}) 

TT X 



AR'''^ix,y;C,m9) = -CF- 

TT X 

R(^Q{x,y;^,7j;g) = ^CF- 

TT X 



y _x 
y _^ 



x/^ < y/v < 1) + 



) + < 








fy- 




[v- 





+ 6{ri-y)e{0<x<^) 
+ d{fi-y)0{O<x<o}, 



(A9) 



(AlO) 
(All) 

(A12) 
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ARQ'^ix, V- 77; g) = ^Nf i 1 1 [5{x/i - y/r^) 9{y < rj) + S{x/^ - y/fj) 9{y > rj) 

-9{Q <x/£,< mm{y/ 77, y/ f)}) 
y; v; 9) = Ai?Q^(:., y- f , r;; .9) + 4 ^iV^ 1 r;^ f ^ + ^ - 2 



-■/7 Wl - + - - ^7 I ^'(0 < < miii{y/77,2//77}). 



(A13) 
(A14) 



To find a formal solution of the evolution equations for double distributions, we proposed in Refs. to 
combine the standard methods used to solve the evolution equations for parton densities and distribution ampli- 
tudes. Hence, let us start with taking the moments with respect to x. Utilizing the property R'^''{x,y; ^,ri; g) = 
i?'^^(x/e,y;l,r?;5)/Cweget 



(A15) 



where F^(y \ fi) is the nth x-moment of F°'{x, y \ fj.) 



(A16) 



The kernels R'^{y, rj; g) and analogous kernels AR'^{y, rj\ g) governing the evolution of G^(?/ 1 /i) are given by 

n+l 



y,r^-g)^AR2Q{y,^-g)^^CF\(y- 



1 



1 



+ 1^ 



n+l 



1 



1 



n+l y ~ 1] 



1 



n + l 1] — y 
1 



e{y>v)-,:Siy-v)-26iy-v) 



9{y < r,) 



■ dz 



1~ z 



(A17) 



TT 



y \ / 4 



n+l Tj ~ y 



e{y < r,) 



y-^y 

7/ ^ ?? 



5{y ~ 1]) 



J3o_ 

2Nr 



dz 
1 - z 



, (A18) 



R-.-iy,v;9)-AROo^y^,,g)+i2^N'\(yY"(^jy-^ ■ ^y-y 

n n + 1 \ \n J \ n n + 2 



(A19) 



9{y <V)+ (I j Oiy > ,) , 



(A20) 



R<^-iy,V;9)-AR<^-iy,,;g)+A^N,^\ly] f ^ - 3^ ] % < ' '^^^ 



TT TT. + 1 1^ \77 



n n + 2 



rj rj 



(A21) 



AR^Q{y, rj; g) ^^C^-l S{y ~ v) " ^ 
TT n \ n+l 



)n+l / ~ \ 

0{y <v)+ f |j 0{y > v) 



(A22) 



R^'^{y,m9) = -Cf- \siy-v) + ^ 

TT n \ rt + 1 



n+l 



9{y < 77) 



n+l 



9{y > 77) 



(A23) 



From Eqs. (A20), (A22) one can derive the following reduction formulas for the nondiagonal kernels: 



dy 



(A24) 
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n^o dy 



(A25) 



The same relations connect the nondiagonal kernels i?*^'^^ ]^QG -^^jth the exclusive kernels V'-^^{y, rj\ g), V^^{y, 77; g) 
given in ref. To understand their structure, one should realize that constructing the nondiagonal QG and GQ 
kernels, one faces mismatching (zp) factors which in the pure C = 1 case are converted into derivatives with respect 
to y. 

It is straightforward to check that all the kernels R'^'^{y,r]; g) (and ARf^'^{y,ri; g)) have the property 

Rt{y,m9)wn{Ti) = Rniv,y;9)wniy), 

where Wn{y) = (yy)"'^^. Hence, the eigenfunctions of the evolution equations arc orthogonal with the weight 

Wnjy) — {yy)"^^, i-e., they are proportional to the Gegenbauer polynomials C^^^^^'^(y — y), see [§4[0 and Refs. 
p2|-p5[ where the general algorithm was applied to the evolution of flavor-singlet distribution amplitudes. 
Expanding the moment functions F^{y \ jj) over the Gegenbauer polynomials C^^^^^(y — y) 



we get the evolution equation for the expansion coefficients 



(A26) 



(A27) 



where FJ^^ are the eigenvalues of the kernels {y , rj; g) related to the elements 7^'' of the usual flavor-singlet 
anomalous dimension matrix 



nk 



QQ -pQG _ „^V^ y 
hi+k ^ nk " 'ri+fe ' ^ nk 



QG 



GQ 



.^GQ 
, hi+k 



^GG 
nk 



^GG . 

in+k ' 



and similarly for the helicity-sensitive quantities AF^^. Namely, 



-Cf 



1 



N+l 



In 



1 



A7; 



GG 
N 



-2Nr 



[N +1){N + 2) 
1 

iV(7VTT) ^ (A^ + 2)(A^ + 3) 

N+l 

3 



{N + 2){N + S) 



N+l 



N+l 

2 ' 



2 ' 



In = 



iV^ + 3iV -f 4 
(iV + l)(Ar + 2) 
N{N + 3) 



N{N + l)iN + 2){N + 3) 



1 



^^^'^~^^(A^-M)(7V + 2) ' "^^^ "''^(iV + 2)(iV-f 3)- 



(A28) 

(A29) 
(A30) 

(A31) 

(A32) 
(A33) 



Let us consider first two simplified situations. In the quark nonsinglet case, the evolution is governed (in helicity- 
averaged case) by j^^i. alone: 



F^'iy I /i) = {yyr+' ^ A"'=C^^/^(y - y) [log(M/A)] 



l27?_?.//3o 



(A34) 



fc=0 



Since 7^*^ = while all the anomalous dimensions 7^"^ with > 1 are negative, only F^^{y \ /i) survives in the 
asymptotic limit fi —> 00 while all the moments F^^{y \ ji) with n > 1 evolve to zero values. Hence, in the formal 
/i — > 00 limit, we have 



F^^{x, y\fi^ 00) S{x)yy 
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i.e., in each of its variables, the hmiting function F^^{x,y \ fj, oo) acquires the characteristic asymptotic form 
dictated by the nature of the variable: S{x) is specific for the distribution functions 37|, while the yy-form is the 
asymptotic shape for the lowest-twist two-body distribution amplitudes For the nonforward distribution 

of a valence quark q this gives 

r^^'-'\X I A* ex.) = 6N,X{1 - X/O/e , 

where Ng is the number of the valence g-quarks in the hadron. 

Another example is the evolution of the gluon distribution in pure gluodynamics which is governed by 7^^^. 
with (3o = llA^c/3. Note that the lowest local operator in this case corresponds to n = 1. Furthermore, in pure 
gluodynamics, 7p*^ vanishes while 7^*^ < if > 2. This means that in the 11 ^ 00 limit we have 

xF'^{x,y \ ^ 00) ^ 30Six)iyy f 

for the double distribution which results in 

J^f (A I ^ c3o) = 30X2(1 - X/C)VC^ 

for the nonforward distribution. In the formulas above, the total momentum carried by the gluons (in pure 
gluodynamics!) was normalized to unity. 

In QCD, we should take into account the effects due to quark-gluon mixing. Diagonalizing Eq.( A27 ), we obtain 
two multiplicatively renormalizable combinations 

Fnk=Fg+»nkF^k (A35) 

where (omitting the nk indices) 



_L_ (^jGG _ ^QQ ^ ^i^^GG _ ^QQ)2 + ^^GQ^QG^ . (A36) 



Their evolution is governed by the anomalous dimensions 



7± = i (^7^^ -I- ± ^{-fGG _ ^QQ)2 ^ 4^GQ^QG j . (A37) 

In particular, 7^^ = and a^Q — 1 which means that = F^f^ + F^ does not evolve: the total momentum carried 
by the partons is conserved. Another multiplicatively renormalizable combination involving F^q and is 



It vanishes in the /i — > 00 limit, and we have 



Since all the combinations F^j^ with n + k > 2 vanish in the /i — s- cx) limit, we obtain 



xF''{x,y\fi^^)^30^^^S{x){yyf ; xFQ{x,y\ ^ 00) ^ 30 ^^^^ S{x){yyf , (A39) 



F'3(x, y I M - 00) -> -30 -^^^^ 5'{x){yyf . (A40) 



In terms of nonforward distributions this is equivalent to 



^f(A|,^oo)^30^^i^|^(l-fV, (A41) 
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Note that both T^iQ and T^{C) vanish in the ji ^ oo hmit. 
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FIG. 2. a) Integration lines in the (a:,y)-plane giving reduction of double distributions F{x,y;t = 0) to usual parton 
densities f{xi) and f{x2)- b) Symmetry line y = {1 — x)/2 for double distributions. 




20 



X, 1^ X, 1 ^ X-!^ !^ X 1 X X, X 1 ^ 



a) b) c) 

FIG. 4. Relation between double distributions F{x,y) and nonforward parton distributions J^(;{X). a) Integration lines 
for three cases: Xi > X — and X2 < C- c) Comparison of integration lines for the nonforward parton distribution 
jFj(X) and usual parton densities f{X), f{X') (shown in 26) and f{X), /(X2) with X2 = X' (shown in 2c). 




FIG. 5. Nonforward parton distributions JT; {X) for different values of the skewedness (" = 0.05 (thin line), = 0.1 
(dashed line), C = 0.2 (dash-dotted line) and C = 0.4 (full line) in the "valence quark oriented" model specified by Eq.(6.11) 
for a = 0.5. 




FIG. 6. Ratio J^l{X)/Xfg{X) vs. logio(^/C) as obtained from the model given by Eq. (6.16) for o) C = 10"^, 6) 
C = 10"^, c) C = 10"* and d) C, = 10"^ with k = 0.48 (solid lines), k = 0.4 (long-dashed lines) and k = 0.30 (short-dashed 
lines). 
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FIG. 7. Parton picture in terms of o) modified double distributions and b) off-forward parton distributions H{x,^). 
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